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20  ABSTRACT  (Continua  on  ravaraa  alda  It  nacaaaary  and  Idantlly  by  block  numbar) 

This  investigation  was  initiated  to  develop  analytic  procedures  for  estimating  aerodynamic 
coefficient  derivatives  for  missiles.  The  analytic  estimates  will  depend  primarily  on  the  geometrical  con- 
figurations of  the  missiles. 

The  problem  of  determining  the  coefficient  derivatives  becomes  reasonably  tractable  for  thin 
airfoil-body  combinations  with  moderate  finite  aspect  ratios  and  flying  at  sonic  speeds.  Starting  with 
the  equations  for  a perfect  gas,  a linearization  of  them  is  achieved  by  assuming  flow  over  a thin  - — 
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20.  ABSTRACT  (Continued): 

prprofile.  A further  assumption  of  a speed  of  Mach  1 gives  rise  to  slender  body  theory.  The  problem 
is  thus  reduced  to  a potential  boundary  value  problem  in  a cross-flow  plane.  Upon  consideration  of 
the  total  momentum  in  a cross-flow  slab,  it  is  found  that  the  resultant  lateral  force  may  be  expressed 
as  a contour  integral  of  the  velocity  potential.  The  effects  of  missile  angle  of  attack  and  control- 
surface  angle  are  incorporated  by  way  of  Neumann-type  conditions  on  the  boundary  contour. 


For  the  special  case  where  the  missile  cross  section  is  a circle  with  midwing,  there  is  an  analytic 
solution  for  the  potential-flow  problem.  For  the  case  where  there  is  an  arbitrary  missile  cross  section 
a computer  program  has  been  developed  which  addresses  the  problem  using  a source  distribution  ap- 
proach. Results  are  given  for  several  sample  cross  sections.  It  is  shown  how  the  cross-flow  results 
may  be  applied  to  a typical  missile  configuration  to  obtain  the  aerodynamic  coefficient  derivatives. 
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AERODYNAMIC  COEFFICIENT  DERIVATIVES  OF  SONIC 
MISSILES  VIA  SLENDER  BODY  THEORY 


INTRODUCTION 


The  Effectiveness  of  Navy  Electronic  Warfare  Systems  (ENEWS)  Program  is  involved  in 
continuing  development  of  simulation  facilities  at  NRL  to  allow  U S.  Navy  investigation  of  the 
interaction  between  anti-ship  missiles  and  shipboard  defensive  electronic  warfare  (EW)  systems. 
Over  the  past  four  years,  the  ENEWS  Program  has  developed  techniques  for  assessing  evolving 
EW  systems  using  mathematical  models  in  simulation  programs  on  digital  computers. 


The  advent  of  this  new  generation  of  models  is  supplemented  by  the  arrival  of  new 
methods  of  missile  analysis  for  accomplishing  the  intent  of  EW  simulations.  To  achieve  realis- 
tic simulation,  it  is  imperative  that  the  aerodynamic  properties  of  the  threat  missiles  to  be 
modeled  must  be  closely  approximated  using  valid  analytic  results.  This  report  forms  a link  in 
the  chain  of  continuing  investigations  intended  to  assess  missile  threats,  current  and  projected. 
This  continuing  assessment  is  an  essential  input  to  the  evaluation  and  enhancement  of  elec- 
tronic warfare  systems  intended  to  counter  missile  threats. 


At  its  current  stage,  this  assessment  program  is  developing  expeditious  and  versatile  pro- 
cedures for  estimating  transonic  aerodynamic  coefficient  derivatives.  These  values  arc  required 
in  computer  simulations  that  incorporate  missile  flight  behavior.  This  information  is  also  useful 
in  determining  outer  limits  for  the  maneuvering  capability  of  a missile,  thus  providing  inputs  to 
bounding  the  potential  capabilities  of  threats.  An  analytic  ability  to  determine  aerodynamic 
coefficient  derivatives  is  essential  for  the  comprehensive  study  of  missile-EW  engagements. 
These  coefficients  depend  upon  the  geometric  configuration  of  the  missile,  including  contro1 
surface  positions,  the  velocity  of  the  air  stream  relative  to  the  missile,  and  the  orientation  o 
the  missile  relative  to  the  air  stream.  The  values  of  the  various  aerodynamic  coefficient  deriva- 
tives may  be  estimated  using  wind-tunnel  tests  on  missile  models,  or  alternatively,  by  means  of 
analytic  calculation.  Although  wind-tunnel  tests  are  preferred  from  the  point  of  view  of  accu- 
racy and  validity,  the  great  expense  involved  in  performing  them  must  be  taken  into  account. 
Wind-tunnel  tests  will  continue  to  be  used  to  determine  the  aerodynamic  coefficient  derivatives 
for  a few  of  the  more  important  missiles.  For  the  multitude  of  others,  it  would  be  preferable  to 
have  and  utilize  a capability  for  making  analytic  estimates  of  these  coefficients. 


This  problem  becomes  tractable  for  missiles  with  moderate  aspect  ratios  at  sonic  speeds. 
In  such  cases  assumptions  may  be  introduced  whieh  greatly  simplify  the  mathematical  descrip- 
tion of  the  problem.  By  treating  a missile  of  moderate  aspect  ratio  as  a thin  profile  and  regard- 
ing the  flow  about  the  missile  as  a perturbation  of  a uniform  flow,  it  follows  that  the  basic  equa- 
tion of  a perfect  gas  can  be  rendered  linear  Requiring,  in  addition,  that  the  flow  be  sonic 
reduces  the  formulaton  to  that  of  slender  body  theory. 
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In  addition  in  satisfying  a number  of  requirements  lot  sonic  missile  analysis.  these  evalua 
lions  would  he  useful  in  missile  analysis  outside  the  some  range  Slemlei  body  theory  may  he 
used  10  Unit  wimr  body  intetference  lactois  that  can  he  used  lo  ileietmme  valid  esiimales  ill  the 
suhsomc  ami  su|*ersomc  values 

Caution  should  he  exercised  however  in  applying  the  results  ol  slender  hod)  iheoiv  1 lie 
direct  results  of  this  iheotv  should  he  tempered  with  some  understanding  of  transonic  flow  and 
its  main  complexities,  otherwise  the  uninformed  usei  could  potentially  he  led  into  attributing 
properties  to  transonic  vehicles  which  are  erroneous 

1 he  theoretical  development  on  which  the  estimation  of  these  coefficient  derivatives  is 
based  will  he  outlined  herein  fins  theory  reduces  the  problem  to  one  ol  solving  a two- 
dimensional  potential  cross  (low  problem  In  the  special  case  where  a cross  section  ol  the  mis 
sile  is  a circle  with  mid  wing,  there  are  analytic  solutions  for  the  potential  (low  problem  l ot 
other  cross  sections,  a computet  program  has  been  developed  m the  course  ol  this  investigation 
that  can  he  used  to  solve  the  two  dimensional  potential  cross  llow  problem 

SIMMARY  OK  I'RIXUU  HK  KOK  COKKK Il'lKNT 
IH  RIV.vriVK  IfKTKRMIN  ATION 

This  report  embodies  a procedure  for  determining  the  sonic  aerodynamic  coefficient 
derivatives  of  a missile,  given  its  geometric  configuration  This  procedure  ts  summarized  and 
certain  equations  are  referred  to  that  will  he  discussed  or  derived  later  in  the  report  According 
to  this  procedure,  certain  cross  section  profiles  are  chosen  at  various  locations  along  the  length 
of  the  missile  Fach  individual  cross-section  profile  is  subjected  to  analysis  using  a computer 
program  designed  to  determine  cross  flows.  To  obtain  input  for  this  computer  program,  it  is 
necessary  to  express  the  cross-section  profile  quantitatively  To  accomplish  this,  the  profile  may 
he  broken  down  into  line  segments  and  arcs,  representing  each  with  a tew  numerical  parame- 
ters Another  input  to  the  computer  program  is  a set  of  boundary  conditions  that  is  determined 
by  body  angle  and  control  surface  deflections  From  these  inputs,  the  cross-flow  program  com- 
putes the  apparent  area,  that  is  a quantity  representing  the  overall  aerodynamic  effect  of  a cross 
section  The  cross-flow  program  represents  the  nuyjor  portion  of  the  computational  effort  in 
determining  the  aerodynamic  coefficients  By  use  of  Fqs.  081-1411.  the  apparent  areas  are 
translated  into  the  aerodynamic  coefficient  derivatives  characteristic  of  the  individual  cross  sec 
tions 

The  acrodvnamtc  coefficient  derivatives  characteristic  ol  the  overall  missile  are  obtained 
from  the  coefficient  derivatives  characteristic  of  cross  sections  using  l-'qs  ) 4 ? > tSH  Among 
the  othet  inputs  to  these  equations  are  the  positions  of  chosen  cross  sections  and  the  moment 
reference  center 

AKRODYNAMIC  < OKKKU  IKN  I DKRIVATIVKS  DKKINKI) 

A number  of  the  more  important  aerodynamic  coefficient  derivatives  will  he  defined  aftei 
the  introduction  of  notations  and  figures  that  describe  and  illustrate  the  mathematical  expres 
sums  and  sign  conventions  used  All  angles  are  expressed  in  degrees 

I he  notations  used  in  defining  the  aerodynamic  coefficient  derivatives  are  presented  next 
and  in  Figs  | and  2 
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(reference  diameter) 
(normal  force) 

(side  force) 

(pitching  torque) 

(yawing  torque) 

(dynamic  pressure) 

(reference  area) 

(free  stream  velocity  ) 
(angle  of  attack) 

(sideslip  angle) 

(elevator  angle) 

(rudder  angle) 

(angular  velocity  of  pitch) 
(density  of  air) 

(yawing  angular  velocity) 


The  pitch  aerodynamic  coefficient  derivatives  considered  are  defined  as  follows: 
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These  sets  of  coefficient  derivatives  may  be  used  to  express  the  normal  force  due  to  angle 
of  attack  and  elevator  angle.  Thus, 

F„  - qS(Cs„a  + Cy,fi).  (6) 


The  equation  for  pitching  torque  due  to  angle  of  attack,  elevator  angle,  and  angular  velo- 
city of  pitch  is 

M - qSb(Cmna  + Cm,br)  + Cm„  9.  (7) 


The  yaw  aerodynamic  coefficient  derivatives  used  herein  are  defined  as  follows: 
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The  yaw  coefficient  derivatives  may  be  used  to  express  the  side  force  due  to  sideslip  angle 
and  rudder  angle  Thus, 

F,  - + C,„5>).  U3> 

Yawing  torque  due  to  sideslip  angle,  rudder  angle,  and  angular  velocity  of  yaw  is 
expressed  as 


N - QSb(C\^  + C.,8,)  + C.> 


U4> 


SLENDER  BODY  THEORY 


In  setting  forth  the  rationale  underlying  the  use  of  slender  body  theory,  it  is  assumed  that 
we  are  dealing  with  the  flow  of  a perfect  gas  about  a solid  body.  The  velocity  potential  4*  for 
the  flow  of  this  compressible  fluid  is  governed  by  the  following  partial  differential  equation  111 
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where  the  spatial  coordinates  are  denoted  by  y.  and  r The  local  velocity  of  sound  c is  given 
by 

. , (U*IJ  U*IJ  U*l:l 

(16) 
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where  the  following  notation  is  used: 


‘o 

(stagnation  value  for  the  speed  of  sound) 

* - eye. 

(ratio  of  the  specific  heats) 

o 

(specific  heat  at  constant  pressure) 

c, 

(specific  heat  at  constant  volume) 
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Next,  the  previously  given  partial  differential  equation  is  linearized.  To  do  this,  consider  a flow 
about  a thin  profile  which  is  parallel  to  the  free  stream  This  flow  may  be  represented  as  a uni- 
form flow  plus  a small  perturbation  In  the  following  equation  the  velocity  potential  <?'  is 
represented  as  a uniform  flow  potential  )'».v  plus  a velocity  perturbation  potential  </> 

<t>  - + <b  (17) 


Substituting  Eq  (17)  into  Eqs.  (15)  and  (16)  results  in  a differential  equation  for  the  velocity 
perturbation  potential  Upon  linearization,  the  equation  becomes 


2 b:<t>  , d:<t>  b:<!> 

a.v2  a.v  b: 2 

(18) 

where 

y2  - 1 - Ml. 

(19) 

VJc„. 

(20) 

and 

•»  i k 1 , . •» 

ci  “ c<5  2 » i- 

(21) 

The  development  that  results  when  the  first  term  in  Eq  (18) 
"slender  body  theory"  This  theory  yields  the  following  equation: 

is  neglected  is  referred  to  as 

b:<t>  b:<b  _ 

a.v*  6;‘ 

(22) 

According  to  this  approximation,  only  flow  perturbations  confined  to  planes  perpendicular  to 
the  body  are  considered  significant 


The  primary  application  of  slender  body  theory  is,  as  the  name  suggests,  to  slender 
bodies  These  are  configurations  whose  lateral  dimensions  such  as  span  and  thickness  are  small 
compared  to  their  length  And,  for  this  type  of  missile,  the  expression  in  the  first 

term  of  Eq.  (18)  may  be  neglected,  giving  rise  to  Eq.  (22). 

Slender  body  theory  may  also  apply  to  sonic  flows  about  non-slender  bodies.  The  small 
disturbance  assumption  is  still  valid  for  a thin  wing-body  combination  with  a moderate  finite 
aspect  ratio,  and  for  sonic  flows  y in  Eq.  18  may  be  neglected,  again  resulting  in  Eq.  (22). 
Under  these  conditions,  slender  body  theory  will  provide  a first-order  approximation  to  sonic 
lifting  flows  [2] . 

Next,  it  is  shown  just  how  estimates  of  the  aerodynamic  coefficient  derivatives  may  be 
obtained  using  slender  body  theory.  Consider  an  infinite  slab  of  air  of  thickness  dx  perpendicu- 
lar to  and  moving  with  the  free  stream  According  to  Eq.  (22),  a two-dimensional  potential 
flow  within  the  slab  is  assumed.  The  problem  of  the  missile  cutting  through  the  slab,  as  illus- 
trated in  Fig.  3,  may  be  treated  as  a plane  flow  with  an  immersed  two-dimensional  region  The 
moving  boundary  contour  of  this  region  provides  boundary  conditions  for  the  potential  flow 
The  effect  of  missile  angle  of  attack  is  to  impart  a uniform  velocity  to  the  contour  within  the 
plane,  as  illustrated  in  Fig  4. 
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Kig  3 — Missile  culling  through  a slab  of  air 
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Fig  4 — Effect  of  body  angle 
on  boundary  conditions 


The  effect  of  control  surface  deflection  is  to  impart  a velocity  to  that  part  of  the  contour 
corresponding  to  the  control  surface  as  illustrated  in  Fig  5.  An  increase  in  size  and  extent  of 
the  missile  cross  section  as  it  cuts  through  the  slab  gives  a point  on  the  contour  an  outward 
velocity  as  shown  in  Fig.  6 


An  additional  boundary  condition  is  provided  by  a zero  fluid  velocity  at  infinity  in  the 
cross-flow  plane.  These  conditions  determine  (except  for  an  additive  constant)  ''-e  two- 
dimensional  velocity  potential  </>  within  the  slab 
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Fig  5 — Effect  of  control  surface  angle 
on  boundary  conditions 
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Fig  6 — Effect  of  expanding  cross  section 
on  boundary  conditions 


The  total  momentum  vector  of  the  flow  d P is  expressed  in  the  following  equation  as  a 
surface  integral  of  the  momentum  density,  over  the  region  exterior  to  the  boundary  contour. 

d P - dx  IS  p V dA,  (23) 

where  the  following  notation  is  used: 

d P (total  momentum  vector  in  slab) 

dx  (thickness  of  slab) 

V (local  velocity  vector  of  fluid  within  plane  slab) 

dA  (element  of  surface  area) 


The  surface  integral  in  Eq.  (23)  may  be  expressed  as  a line  integral  around  the  boundary  con- 
tour as  follows: 

d?  - dx  p $ <t>  di.  (24) 

where  di  is  an  outward  normal  vector  representing  an  element  of  arc  length  along  the  boundary 
contour. 
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The  force  vector  r/F  on  a missile  cross-sectional  element  will  be  equal  and  opposite  to  the 
rale  of  increase  of  momentum  in  the  slab. 


This  is  expressed  as 


r/F  - — — (d P). 
dt 


(25) 


Taking  into  consideration  that  the  missile  is  cutting  through  the  slab  with  a velocity  Eq. 
(25)  becomes 


d F - - -7-  (</ P). 

dx 


(26) 


Slender  body  theory  has  greatest  validity  for  wings  with  monotonically  increasing  span  [3). 
In  calculating  particular  aerodynamic  coefficient  derivatives,  the  analyst  will  riot,  in  most  cases, 
be  dealing  with  the  simple  case  of  a gradually  expanding  contour.  The  body  may  be  tapered  so 
that  it  gives  a contracting  contour,  or,  as  at  the  trailing  edge  of  a wing,  there  may  be  a discon- 
tinuous reduction  in  the  extent  of  the  contour. 

In  order  to  deal  with  these  situations,  it  will  be  taken  as  a hypothesis  that  Eq.  (26)  gives 
the  r/F  only  when  the  magnitude  of  d P is  increasing  along  the  axis  of  the  missile.  Otherwise, 
r/F  will  be  taken  to  be  zero. 


To  obtain  the  total  lateral  force  vector  F on  a missile,  the  force  vector  dF  on  a cross  sec- 
tional element  is  integrated  along  the  length  of  the  missile,  the  coordinate  x being  used  to 
denote  positions  along  that  length  down  stream  from  the  nose.  Consider  an  interval  along  the 
length  of  the  missile  over  which  the  magnitude  of  dP  is  increasing.  The  total  lateral  force  on 
this  interval  is  expressed  as 

\-\b 

F-  f r/F.  (27) 

'-'4 

Here,  x4  and  xH  are  the  end  points  of  the  interval.  The  integration  is  carried  out  after  substi- 
tuting Eq.  (26)  into  Eq.  (27),  thus  obtaining 


Substituting  Eq.  (24)  into  Eq.  (28)  results  in  the  following  equation  for  the  total  lateral  force 
on  the  interval: 


F 


- p <t>  di 


(29) 


A particular  missile  cross  section  with  a specified  direction  of  motion  within  a cross-flow 
slab  could  be  characterized  using  the  conventional  quantity  apparent  mass.  However,  we  shall 
define  a related  quantity,  the  apparent  area  of  the  cross  section,  a more  convenient  unit  to 
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express  this  concept.  We  define  this  vector  quantity  according  to  the  following  equation,  where 
the  surface  integral  is  over  the  region  exterior  to  the  boundary  contour; 


A - 


iffy* 


(30) 


The  expression  Vc  is  a velocity  naturally  associated  with  the  motion  of  the  contour  in  the 
cross-flow  plane.  This  will  generally  be  either  the  uniform  velocity  of  the  contour  due  to  body 
angle  or  the  velocity  of  that  portion  of  the  boundary  representing  a control-surface  deflection. 
It  is  given  as  follows. 


V,  - 


180 


(31) 


Here,  rj  is  the  body  angle  or  control-surface  deflection.  The  velocity  K,  may  also  be  due  to  the 
angular  velocity  of  pitch  or  yaw  denoted  by  <u.  In  this  case  it  is  given  by 


Vc 


ir 

180 


\x  - xf|  w. 


(32) 


Here,  x(  denotes  the  position  of  the  moment  reference  center  measured  downstream  from  the 
nose  and  chosen  to  coincide  with  the  center  of  gravity. 

The  significance  of  the  apparent  area  is  that  the  apparent  mass  of  the  cross  section  is 
given  by  p|A|</x,  and  the  momentum  of  the  cross  flow  is  given  by  pVcKdx.  The  apparent  area 
may  be  expressed  as  the  following  contour  integral; 

A“  (33) 

r f 

Note  that  the  forgoing  expression  does  not  depend  on  the  magnitude  of  Vc  since  the  boundary 
condition  values  and,  thus,  the  contour  integral  are  proportional  to  this  velocity. 


The  transverse  force  and  the  moment  of  force  on  a missile  may  be  expressed  in  terms  of 
the  apparent  areas  of  the  missile  cross  sections.  The  transverse  force  due  to  body  angle  or 
control-surface  deflection  on  an  interval  taken  along  the  length  of  a missile  where  the  apparent 
area  is  increasing  is: 


F 


(34) 


The  moment  of  force  on  this  interval  is  given  by 

*-*• 

M - / (x  - xc)  IxdF.  (35) 

y-x* 

Here,  I is  a unit  vector  pointing  downstream. 

Substituting  Eq.  (34)  in  Eq.  (35)  yields  the  following  for  the  moment  of  force  on  this  interval: 


10 


NRL  REPORT  8203 


p Vi  tj  f (V  - O lx</A.  (36) 

In  a similar  manner,  Eq.  (37)  is  obtained  for  the  moment  of  force  due  to  angular  velocity  of 
pitch  or  yaw: 

(if.i.  f (v  - v,)1  IxdA.  (37) 

APPLICATION  TO  A TYPICAL  MISSILE  CONFIGURATION 

A method  of  estimating  the  aerodynamic  coefficient  derivatives  of  a typical  missile 
configuration  is  described  here.  Top  and  s*de  views  of  the  missile  configuration  under  con- 
sideration are  shown  in  Figs.  7 and  8 Various  points  along  the  length  of  the  missile  have  been 
numbered  as  follows: 


Fi|  8 — Side  view  of  typical  missile  configuration 
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1 (a  point  where  the  cross  section  of  the  nose  attains  a maximum), 

2 (beginning  of  wing), 

C (moment  reference  center), 

3 (trailing  edge  of  wing), 

4 (beginning  of  tail), 

5 (trailing  edge  of  elevator  and  rudder). 


The  value  of  v at  the  / th  point  along  the  length  of  the  missile  is  denoted  by  x,. 

One  aspect  of  carrying  out  the  method  for  determining  the  aerodynamic  coefficient  deriva- 
tives is  to  determine  apparent  areas  for  the  cross  sections  at  a number  of  points  along  the 
length  of  the  missile  Notation  is  given  below  for  apparent  areas  as  given  by  Eg.  (33)  for  the 
cross  section  located  at  the  station  denoted  by  v.  For  angle  of  attack  the  apparent  area  is 
denoted  by  Av„(v),  where  in  this  case  17  in  F.q  (31)  represents  the  angle  of  attack  of  the  mis- 
sile, a.  The  effect  of  this  a for  small  angles  is  to  give  the  boundary  contour  a rigid  motion 
with  velocity  F.,  as  illustrated  in  Fig  4,  thus  determining  the  boundary  conditions  for  <i> 

For  elevator  deflection  the  apparent  area  is  denoted  by  A v*(x),  where,  in  this  case,  17  in 
Eq.  31  represents  the  elevator  deflection  6,..  The  effect  of  this  for  small  deflections  is  to  contri- 
bute to  that  part  of  the  boundary  representing  the  control  surfaces  a velocity  Fc,  while  the 
remainder  of  the  contour  assumes  a zero  velocity.  This  is  illustrated  in  Fig  S and  provides  the 
boundary  conditions  for  <t>. 

For  side  slip  angles  the  apparent  area  is  denoted  by  A )(j(.v),  where,  in  this  case,  17  in  Eq. 
31  represents  the  angle  of  sideslip  of  the  missile  The  effect  of  this  0 for  small  angles  is  to 
give  the  boundary  contour  a rigid  motion  with  velocity  V.s  thus  determining  the  boundary  con- 
ditions for  <t> 

For  rudder  deflection  the  apparent  area  is  denoted  by  Af#(.v),  where,  in  this  case,  17  in 
Eq.  (31)  represents  the  rudder  deflection  8 >•.  The  effect  of  this  for  small  deflections  is  to  con- 
tribute to  that  part  of  the  boundary  representing  the  control  surface  a velocity  Fr,  while  the 
remainder  of  the  contour  assumes  a zero  velocity.  This  provides  the  boundary  conditions  for 


The  following  equations  introduce  cross  section  coefficient  derivatives.  These  are  related 
in  a simple  way  to  the  apparent  areas  and  are  useful  in  expressing  the  aerodynamic  coefficient 
derivatives  of  the  missile.  Thus: 


Cv„  <*)  “ - 

9 h 

(38) 

” — 

*».(*)!• 

(39) 

Cr»  (*)  - ~ 

WsA»(xhk- 

(40) 

Cy»(x)  -- 

90  SAr,ix)  k 

(41) 

and 
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Here,  j and  k are  unit  vectors  in  the  cross-flow  plane  and  in  the  directions  of  FN  in  Fig.  1 and 
Fy  in  Fig.  2,  respectively. 


Formulas  for  estimating  the  aerodynamic  coefficient  derivatives  of  the  missile 
configuration  are  listed  next. 


CNa  “ N n ^Vq  'X4)  C ,V a (x5) 

Cvs  “ Cv»(x5) 


C — 

b 2 


- mb 


x\  xi  XS 

J + / + / | (x  - xc)  dC\a  (x) 

0 Jtj  x4 

*'  *1 

f + J + / (x  - xf)2  dC\a(x) 

0 Xj  Xt  J 

- ~ -r  jf  (x  - x(.)  dCNh(x) 


Cyp-  C yp  (x  3)  — Cyp(x  4)  + Cyp(x  5) 
Crs  ” C (x  5) 

f* 

1 


C^“-T 


c - — — 
*2 


J + J + / (x  " XJ  dCYfiM 

0 x j x4j 

X,  Xj  Xjj 

/+/+/(*-  X(.)  2 dCyp  (x) 
0 Xj  X4J 

1 f 

C"*  “ ~ ^ J (jr  " X<)  4/(7 r# 


(42) 

(43) 

(44) 

(45) 

(46) 

(47) 

(48) 

(49) 

(50) 

(51) 


In  obtaining  the  normal  force  due  to  body  angle,  two  intervals  are  considered.  One  of 
these  extends  from  the  nose  tip  to  the  trailing  edge  of  the  wing  and  the  other  extends  from  the 
beginning  of  the  tail  to  the  trailing  edge  of  the  elevator.  Using  Eqs.  (1),  (34),  and  (38),  the 
coefficient  derivatives  for  the  normal  force  due  to  body  angle  are  expressed  by  Eq.  (42).  In  a 
similar  manner  the  coefficient  for  side  force  due  to  side-slip  angle  is  expressed  by  Eq.  (47). 

The  normal  force  due  to  elevator  deflection  is  obtained  as  the  force  on  the  interval 
extending  over  this  control  surface.  Using  Eqs.  (2),  (34),  and  (39),  the  coefficient  derivative 
for  this  force  is  expressed  by  Eq.  (43).  Similarly,  the  coefficient  derivative  for  side  force  due  to 
rudder  deflection  is  expressed  by  Eq.  (48). 

In  obtaining  the  pitching  moments  due  to  body  angle  and  angular  velocity,  three  intervals 
are  considered.  The  first  extends  over  the  nose,  the  second  over  the  wings,  and  the  third  over 
the  horizontal  stabilizer.  Using  Eqs.  (3),  (36),  and  (38),  the  coefficient  derivative  for  the 
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pitching  moment  due  to  body  angle  is  expressed  by  Eq.  (44).  Using  Eqs.  (4),  (37),  and  (38), 
the  coefficient  derivative  Tor  the  pitching  moment  due  to  angular  velocity  of  pitch  is  expressed 
by  Eq.  (45). 

In  a similar  manner  the  coefficient  for  yawing  moment  due  to  side-slip  angle  is  expressed 
by  Eq.  (49),  and  the  coefficient  for  yawing  moment  due  to  angular  velocity  of  yaw  is  expressed 
by  Eq.  (50). 

The  pitching  moment  due  to  elevator  deflection  is  obtained  as  the  moment  on  the  interval 
extending  over  this  control  surface.  Using  Eqs.  (5),  (36),  and  (39),  the  coefficient  derivative 
for  this  force  is  expressed  by  Eq.  (46).  Similarly,  the  coefficient  for  yawing  moment  due  to 
rudder  deflection  is  expressed  by  Eq.  (51). 

It  is  apparent  from  Eqs.  (42),  (43),  (47),  and  (48)  that  the  coefficient  derivatives  for  the 
normal  and  side  forces,  namely,  CSa,  C*#,  C>()  and  C>,,  are  dependent  on  the  cross  flows  at 
the  three  stations  along  the  length  of  the  missile  numbered  3,  4,  and  5.  The  other  coefficients 
involve  integrations  and  are  evaluated  as  approximating  sums.  This  evaluation  requires  analyz- 
ing a number  of  cross  flows  taken  along  the  intervals  of  integration. 

THE  CIRCLE-W1TH-MIDWING  CROSS  SECTION 

In  the  special  case,  where  the  cross-section  contour  of  a missile  of  angle  of  attack  a is  a 
circle-with-midwing,  the  potential  cross-flow  problem  has  an  analytic  solution.  The  boundary 
condition  for  a problem  of  this  type  is  illustrated  in  Fig.  9.  The  following  notation  is  also  indi- 
cated: 

R ( x ,)  (radius  of  circular  cross  section  of  missile) 

s(x,)  (semi-span  of  wing) 
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For  this  case,  the  cross-section  coefficients  (4|  are  given  by 


RHx)  + 


R\x)  | 
s’(.v)  ) 


(52) 


Applying  this  formula,  the  pitch  coefficient  CVo  may  be  estimated  for  a missile  with  the 
configuration  shown  in  Figs  7 and  8,  and  which  has  circle-with-midwing  cross  sections.  This 
coefficient  derivative  is  given  by  the  formula: 


f v«.  “ 


sJ(.v,)  - R2(x j)  + 


90  S 

- R :(.v4)  + - RHxs)  + 


R\x>) 
sJ(x,) 
R\xd 


S2  Uj) 


(53) 


The  formula  given  by  Eg  (52)  may  also  be  applied  in  determining  the  coefficient  derivatives 
and  C’m<r  This  is  accomplished  by  substituting  into  Eqs  (44)  and  (45). 

A COMPUTER  PROGRAM  FOR  ARBITRARY 
CROSS  SECTIONS 

A computer  program  has  been  developed  which  solves  the  potential  flow  problem  for 
shapes  of  the  type  encountered  with  missile  cross  sections.  This  program  also  determines  the 
apparent  area  of  the  cross  section  through  evaluation  of  the  integral  expressed  in  Eq.  (33). 

A missile  cross  section  can  generally  be  represented  by  a region  representing  the  body 
along  with  line  segments  representing  wings  and  various  other  airfoils  The  boundary  of  this 
configuration  will  be  the  internal  boundary  for  the  potential-flow  problem  Body  angle  or 
control-surface  deflection  is  represented  by  specifying  the  normal  gradient  along  this  boundary 
It  is  further  assumed  that  the  gradient  of  the  potential  goes  to  zero  at  infinity 

The  calculation  is  based  upon  source  distribution  methods  similar  to  those  which  have 
been  used  by  J.  L Hess  and  others  |5  — 7)  in  several  computer  programs  With  these  methods, 
a general  potential  function  is  represented  as  a linear  sum  of  a fundamental  system  of  potential 
functions,  each  of  which  is  derived  from  a source  of  some  type.  The  source  may  be  a point 
monopole  or  dipole  or  a continuous  distribution  of  one  of  these  types. 

Thus  we  may  represent  a general  potential  function  as 

\ 

<M\|.  Vj)  - £ <t>  (,V|,  vd.  (54) 

-i 

Here,  x , and  x • are  coordinates  in  the  cross-flow  plane  The  expression 

<t>  (.V|.  ,Vj)  i — l N 

represents  the  fundamental  system  of  source-derived  functions.  Each  of  these  functions  takes 
on  a weight  denoted  by 

C, 


i - 1 N. 


rih-Hr-*; 


▼ 
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We  may  derive  the  velocity  field  from  the  potential  by  taking  the  gradient  of  the  latter,  as 
shown  by  the  following  equations: 


and 


V(x„  x2) 


d<t>  d<t> 
fix,  ' dXj 


(55) 


V,(x„  xj) 


fi<6,  d<t>, 

fix  | fix  2 


(56) 


Applying  these  equations  to  Eq.  (54)  results  in  the  following  general  expression  for  the  velocity 
field: 


V(x,.  x2)  - £ C,  V,(x„xj). 

-i 


(57) 


A set  of  N points  are  chosen  on  the  boundary  of  the  contour  as  locations  for  imposing  the 
boundary  condition.  The  coordinates  of  these  points  are  denoted  by  x, J and  x2 •'  where 

(j  - 1 S).  Also,  the  unit  normal  to  the  boundary  at  each  of  these  locations  is  denoted  by 

n,  where  ( j - 1 N).  As  a boundary  condition,  the  normal  velocity  W,  is  specified  at  each  of 

these  points.  The  imposition  of  this  condition  is  expressed  by 

ZC,  V,(x,>  xj4  n,.  (58) 

>-! 

This  relation  represents  A simultaneous  linear  equations  for  the  weight  factors  C,.  Solving  this 
system  of  equations  determines  the  value  of  each  C„  whereupon  the  velocity  potential  <t>  and 
the  velocity  V may  be  calculated  from  Eqs.  (54)  and  (57).  The  apparent  area  vector  may  be 
determined  by  evaluating  Eq.  (33)  with  a numerical  integration  process.  The  general  represen- 
tation of  the  velocity  potential  is  based  on  its  expression  in  terms  of  sources  This  may  be  done 
in  a variety  of  ways.  Many  investigators  (See  Ref.  6 for  citation  of  additional  refs.)  have,  for 
instance,  utilized  continuous  density  distributions  on  the  body  surface. 


For  this  investigation,  the  body  has  been  treated  differently  from  the  airfoils.  For  the 
body,  point  sources  located  within  and  at  some  distance  from  the  boundary  contour  have  been 
utilized.  For  a point  source,  the  velocity  potential  and  the  velocity  at  the  point  (x,.  x2)  are 
given  by  the  formulas 

<t>  - ~ In  [( v , - .v,)2  + (x2  - ,v2)2J  (59) 

and 


(X|  — v,.  -V2  ~ v ?) 
(*i  ~ + (v2  ~ >' 2*2 


(60) 


Here,  v,  and  are  the  coordinates  of  the  point  source.  For  airfoils  which  are  represented  as 
line  segments,  continuous  dipole  distributions  are  utilized.  The  line  segment  representing  a 
wing  is  divided  into  small  elemental  segments,  and,  on  each  of  these,  the  dipole  distribution  is 
represented  by  a quadratic  function.  For  each  elemental  segment,  we  denote  its  length  by  As 
and  introduce  coordinates  x,  and  Xj  such  that  the  segment  lies  on  the  v,  axis  centered  about 
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(he  origin.  To  accomplish  the  general  quadratic  representation,  three  fundamental  potential 
functions  are  introduced  corresponding  to  each  elemental  segment.  These  three  functions  are 
derived  from  dipole  distributions,  the  first  of  which  is  constant  while  the  second  and  third  are 
proportional  to  x,  and  x,J,  respectively. 

For  a constant  dipole  distribution  on  a line  segment,  the  velocity  potential  and  the  velo- 
city at  the  point  x2)  are  given  by  the  following  formulas: 

As  f 2 


4>c  - Incl  |x,J  + xj:  - — | , As  x2 


(61)’ 


and 


V<“ 


As 

2 

2 

As 

i 

As 

2 

X,+  y 

+ xj 

X|-  T 

+ x{ 

(62) 


For  a dipole  distribution  which  is  proportional  to  X|,  the  velocity  potential  and  the  velocity  are 
given  by 

X|  <t> , - x2  <b 


and 


V,  - 


<t>,  - 


l*c.  - 


As 

+ X,  V.  - XjU 


As 


(63) 


(64) 


where  we  have  introduced  the  following  auxiliary  quantities: 


- 2 


and 


II  - 


X|  + 


As 


As 

x,  + y 

} 

+ *i 

i 

As 

Xl_  T 

+ xj 

*2 

X\  - - 

(65) 


• v. 


, As 

Xl+  T 


+ x i 


X I - 


As 


+ 


(66) 


For  a dipole  distribution  which  is  proportional  to  X|J,  the  velocity  potential  and  the  velocity  are 
given  by 


•The  function  Incl  w is  defined  us  the  angle  of  inclination  of  w in  radians  where  rr  < Incl  w ^ rr 
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Vt>- 


<t>o  “ 


(x2  - a/)  d>  , - 2 x | x 2 <l>  + As  .vj 


As2 


2 _Aj|  *,  + (x,2  - x/>  V e - 2 0 |A  | - 2 x,  xj  U + 


*2 


As 


As2 


(67) 


(68) 


In  order  to  determine  the  N weight  factors  C„  an  equal  number  of  linear  equations  are  esta- 
blished Each  of  these  equations  are  generated  by  specifying  the  normal  velocity  W,  at  a point 
on  the  boundary  Those  points  where  the  normal  velocities  W , are  specified  shall  be  referred 
to  as  control  points.  For  each  monopole  source  chosen  internal  to  the  body,  a control  point  is 
selected  on  the  body  boundary.  For  each  elemental  segment  of  a wing,  three  fundamental 
functions  are  derived  This  requires  the  determination  of  three  corresponding  weights  (' 
Therefore,  on  each  such  elemental  segment,  three  control  points  are  selected 


In  formulating  a particular  calculation,  choices  are  made  of  the  control  points  and  ol  the 
sources  from  which  the  fundamental  system  of  basic  functions  are  derived  No  definite  pro- 
cedure has  been  established  for  making  these  choices  in  a way  which  approaches  the  optimal. 
Rather,  making  satisfactory  choices  requires  the  use  of  good  intuitive  judgment  in  combination 
with  knowledge  learned  from  a moderate  amount  of  computational  experimentation. 

Using  a circle-with-midwing  profile,  a certain  amount  of  computational  experimentation 
was  performed.  For  this  profile,  analytic  answers  are  available  for  comparison.  Figure  10  illus- 
trates how  the  sources  were  chosen  for  this  typical  cross-flow  calculation.  The  conclusions 
drawn  from  this  undertaking  are  assumed  to  have  general  validity  for  bodies  with  wings. 


MONOPOLE 
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One  conclusion  deduced  is  that  in  order  to  obtain  accurate  results  it  is  necessary  to  let  the 
line  segment  representing  an  airfoil  penetrate  a distance  into  the  body,  as  shown  in  Fig.  10. 
When  this  line  segment  was  allowed  to  stop  exactly  at  the  body  boundary,  a good  approxima- 
tion to  the  true  solution  was  unobtainable. 

Another  observation  is  that  a modest  increase  in  accuracy  is  obtained  by  letting  the  size  of 
the  elemental  segments  become  smaller  as  the  tip  of  the  wing  is  approached. 

The  information  defining  a particular  cross-flow  problem  is  summarized  next.  Associated 
with  the  body  is  the  position  of  each  interior  monopole  along  with  information  regarding  ele- 
ments on  the  body  boundary.  The  latter  consists  of  the  location  and  orientation  of  each  ele- 
ment along  with  control-point  positions.  Associated  with  an  airfoil  is  the  positon  and  orienta- 
tion of  the  line-segment  representation  along  with  information  regarding  its  partitioning  into 
elemental  segments,  each  elemental  segment  has  three  control  points  designated  on  it.  Also 
included  among  the  given  information  is  the  normal  velocity  W,  which  is  specified  at  each  con- 
trol point  on  a body  boundary  or  airfoil. 

Rather  than  inputting  the  given  information  with  enumerative  lists,  various  portions  of  it 
are  generated  by  subroutines  using  analytic  formulas.  For  example,  the  program  contains  a 
subroutine  which  generates  all  the  necessary  defining  information  for  an  airfoil.  Another  sub- 
routine does  this  for  a circular  arc.  A number  of  subroutines  such  as  these  can  be  used  to  piece 
together  an  entire  cross-section  contour.  Upon  execution,  the  computer  program  calculates  the 
value  of  (he  velocity  potential  4>  and  the  velocity  V at  each  control  point.  It  also  calculates  the 
weights  C,  of  the  various  fundamental  functions.  In  addition,  it  calculates  the  apparent  area 
vector  A as  given  by  Eq.  (33). 

TEST  CALCULATIONS 

In  order  to  test  the  computer  program  for  accuracy  and  computation  time,  the  program 
was  executed  on  a CDC  6600  computer.  Profiles  were  run  for  which  answers  could  be  obtained 
independently  using  analytic  formulas.  The  profiles  used  were  the  circle,  circle  with  midwing, 
and  wing  only.  These  are  all  special  cases  of  the  general  circle-with-midwing  problem  for  which 
an  analytic  result  (Ref.  4)  is  available  for  comparisor  The  general  expression  for  the  apparent 
area  magnitude  A due  to  a vertical  cross  flow  is  given 

A - 7r(sJ-  *j+  RVs1)  (69) 

Here,  R is  the  radius  of  the  circle  and  s is  the  semi-span  of  tne  wing,  as  indicated  in  Fig.  9. 
Analytic  comparison  values  of  the  apparent  area  for  the  profiles  used  in  the  test  computations 
were  obtained  by  using  the  following  three  sets  of  values  in  Eq.  69.  The  values  used  for  R and 
s were:  R - s - 1 for  the  circle;  R - 1/2,  r - 1 for  the  circle  with  midwing;  and  R - 0,  s - 1 
for  the  wing  only.  For  the  circle-with-midwing  profile,  test  calculations  were  performed  using 
various  numbers  of  fundamental  function  generating  elemental  sources.  These  will  be  referred 
to  as  elements.  The  results  of  the  test  calculations  using  24,  42,  84,  and  168  elements  are  sum- 
marized in  Table  1.  In  each  case,  the  number  of  elements  is  broken  down  into  the  number  of 
monopole  types  representing  point  monopoles  within  the  body,  and  the  number  of  dipole  types 
representing  continuous  dipole  distributions  on  an  elemental  segment  of  a wing.  This  table 
indicates  the  improvement  in  accuracy  as  the  number  of  elements  is  increased.  The  0.27% 
error  for  42  elements,  however,  should  be  regarded  as  unusually  low  due  to  chance  and  not 
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really  fitting  the  trend  Table  1 also  indicates  the  rapid  increase  in  computation  lime  as  the 
number  of  elements  is  increased  It  is  apparent  that  a reasonable  combination  of  accuracy  and 
computation  time  is  obtained  when  the  number  of  source  elements  is  on  the  order  of  84  In 
order  to  test  the  program  for  the  effect  of  varying  the  contour  configuration,  calculations  were 
performed  for  the  circle,  circle-with-midwing.  and  wing  only  using  84  elements  in  each  case 
The  results  are  shown  in  Table  2.  indicating  in  each  case  the  individual  number  of  monopole 
and  dipole  elements  used 


Table  I Circle  With  Mjdwing 


APPARENT  AREA  (ANAlVTlCl  2 5525 


No.  of 
ITenients 

No.  of 
Monopole 
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No.  of 
Dipole 
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(Computed) 

Prior 
('•' ) 
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mm 

mm 

18 

2.3257 

8.80 

0.072 

Kfl 

30 

2.5504 

0.27 

2.5o 

84 

■ ^ i 

<i() 

2.5638 

0.44 

12.5 

168 

48 

120 

2.5533 

0.031 

78.1 
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It  is  observed  in  Table  2 that  near-perfect  accuracy  is  obtained  with  the  circular  profile 
(assumed  typical  of  a body-only  type  of  contour).  The  least  accuracy  or  highest  error  (0  44%), 
is  obtained  using  the  circle-with-midwing  This  profile  (assumed  typical  of  a wing-body  combi 
nation)  is  the  case  most  useful  for  revealing  any  possible  error  It  is  lound  that  an  intermediate 
degree  of  accuracy,  or  next  highest  error  (0  1 1%)  is  obtained  with  the  wing-only  profile 

SAMPLE  CALCULATIONS 

This  computer  program  is  of  general  utility  for  performing  cross-flow  calculations  for 
almost  any  shape  that  would  be  encountered  with  a missile  cross  section.  Sample  calculations 
have  been  performed  for  a few  cases  where  analytic  results  aio  not  available  Apparent  areas 
are  calculated  for  boundary  conditions  representing  pi'ch.  elevator  deflection,  yaw.  and  rudder 
deflection  The  results  of  these  calculations  are  given  in  Tables  3.  4,  and  5.  In  Table  3,  a 
rudder  has  been  added  to  the  circle-with-midwing  configuration  Table  4 shows  results  for  a 
configuration  consisting  of  simply  a circular  body  and  a rudder  Results  for  a tri-tail 
configuration  are  given  in  Table  5. 


Table  3 - Circle  With  Midwing  and  Rudder  (114  Elements) 
Type  of  Control  Configuration  Apparent  Area 


Table  4 


Type  of  Control 


Yaw 


Rudder 

Deflection 


Pitch 
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Circle  With  Rudder  (114 


Configuration 


Elements) 

Apparent  Area 
Magnitude 


1.62 


0.538 


0.785 
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Table  5 — Tri-Tail  Configuration  (114  Elements) 


Type  of  Control  Configuration 


Apparent  Area 
Magnitude 


Pitch 


Elevator 

Deflection 


Yaw 


Rudder 

Deflection 


2.14 


1.00 


2.14 


0.579 
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CONCLUSIONS  AND  IMPLICATIONS  FOR  FURTHKR  EFFORT 

Expeditious  and  versatile  procedures  for  estimating  sonic  aerodynamic  coefficient  deriva- 
tives have  been  developed.  This  endeavor  was  based  upon  the  generation  of  a computer  pro- 
gram to  solve  the  potential  cross-flow  problem  This  program  is  similar  to  certain  other  existing 
ones  but  has  some  necessary  differences  The  development  is  based  on  slender  body  theory 
and  gives  estimates  which  can  be  applied  to  missiles  flying  at  sonic  speeds.  These  estimates  are 
also  helpful  in  determining  approximations  for  the  coefficient  derivatives  in  the  subsonic  and 
supersonic  ranges. 

There  is  a considerable  need  for  these  aerodynamic  coefficient  estimates  as  inputs  to  com- 
puter and  hardware  simulations.  These  simulations  deal  with  a wide  variety  of  analyses  of  the 
expected  flight  behavior  of  actual  missiles  under  a variety  of  conditions.  An  extensive  effort  is 
needed  to  investigate  the  implications  of  certain  postulated  characteristics  for  future  missiles. 
This  effort  would  involve  working  out  a missile  design  within  the  constraints  of  the  postulated 
external  characteristics.  The  resulting  design  would  be  extremely  useful  in  bounding  the  threat 
by  indicating  the  outer  limits  of  a missile’s  maneuvering  capability.  Aerodynamic  coefficients 
based  on  this  design  would  provide  valuable  inputs  for  future  simulations  which  could  be  used 
to  study  the  expected  outcome  of  missile-ship  engagements  involving  newer  missile  types. 

A capability  for  analytic  determination  of  aerodynamic  coefficient  derivatives  is  essential 
for  any  comprehensive  investigation  of  missile-ship  engagements.  It  is  needed  not  only  for  the 
analysis  of  existing  missile  threats  but  also  for  the  examination  of  possible  future  missile 
threats. 
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